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The latest search for massive long-lived hadronising particles with the ATLAS detector is 
presented. The search is conducted with the inner detector and an integrated luminosity 
corresponding to 2.06 fb"^ at = 7 TeV. For Split-SUSY scenarios gluino masses below 810 
GeV have been excluded. 



1 Introduction 



> 

, This paper presents the latest search for massive long-lived hadronising particles with the ATLAS 

'nT ' detector at the The discovery of new massive long-Hved particles at the LHC would be of 

' fundamental significance^. It was the discovery of a stable massive particle that opened the field 

. of particle physics, and for many years electrons, protons, neutrons and muons were the only 

Q I known type of elementary particles, all of which could be detected due to their long-lived nature. 

CN ' With the advent of accelerators the discovery of long-lived hadrons heralded the development 

of quantum chromodynamics, the theory of the strong interaction. In that sense history is a 
primary motivation for the search for new long-lived particles. Furthermore multiple new models 
^ I predicts potentially long-lived particles. The search presented here investigates specific SUSY 

^ ' scenarios that predict massive long-lived coloured particles in the form of colour-octet gluinos 

- - - (g). Due to colour confinement it is expected that the heavy sparton will hadronise with standard 

model partons into a colour neutral object called an i?-hadron^. In i?-parity conserving SUSY 
models this bound state is kept stable. The sparton itself is assumed stable or at least long-lived 
by other mechanisms, for instance requiring a highly suppressed decay channel^^l. 

The primary search strategy for these objects relies on them being massive compared to the 
available collision energy. This in turn leads to the assumption that the particles must move at 
low velocities compared to the speed of light, thus behaving differently than most other directly 
observable particles produced. If the slow moving particles carry electric charge the particle can 
be detected in multiple sub systems in ATLAS In general two complementary methods both 
sensitive to particle velo city ; time of flight and anomalous energy deposition are used. 

In previous searches ATLAS reported results based on a combination of sub-detectors, 
Pixel based specific energy loss dE/dx, Calorimeter based time-of-fiight (ToF) and ToF in 
the Muon spectrometer as well. The analysis presented here only uses the Pixel based dE/dx 
method in order to remain insensitive to hadronic interactions in the calorimeters that can cause 
an exchange of light quarks making the i?-hadron electrically neutral and invisible to further 
detection. 



2 Models and datasets 



The analysis is based on LHC pp collision data at ^/s = 7 TeV, recorded in the period between 
March and August 2011. This corresponds to an to an integrated luminosity of 2.06 fb^^ after 
data quality cuts. Signal samples have been generated in a range of masses from 200 to 1000 GeV 
in 100 GeV intervals. The analysis is insensitive to lower masses due to trigger constraints and 
large backgrounds. The R-Hadron samples are pair produced gg simulated with Pythia The 
initial hadronisation of gluinos into R-Hadrons is done with a custom version of the Lund string 
hadronisation model. The particles are always produced stable in this scenario. Simulating the 
R-Hadron interaction with material in the detector is handled by dedicated Geant4 routineJ^^. 
Because of the limited knowledge regarding the possible hadronic states, three different models 
are used to cover various model assumptions^^. In all Monte Carlo samples, the primary collision 
event is overlaid with minimum bias simulated events to model the pile-up conditions in data. 



3 Pixel-based mass estimation 



The silicon pixel detector is the part of the tracker situated closest to the beam pipe in ATLAS 
and covers the radial distance 50.5 < r < 122.5 mm in the pseudorapidity region |ry| < 2.5. The 
detector consists of three layers parallel to the beam pipe in the barrel region and perpendicular 
in in the end-cap regions. In total the pixel tracker contains 80 million pixels. The pixel tracker 
is capable of reading out a "time over threshold" (ToT) value that corresponds to the amount 
of charge collected in each pixel. The threshold is set at a minimum value of 3200 ±80 e~ . The 
maximum charge collected per pixel is limited by the dynamic range of the read out electronics 
(8 bit ADC) to roughly 8.5 times the energy deposited by a minimum ionising particle (MIP). 
Depositions higher than that causes a saturation leading to a mismeasurement. The maximum 
value limits the sensitively to slow particles with velocities higher than f3 > 0.3. To arrive at 
an estimate of the specific energy loss dE/dx the charge deposition is measured in all pixels in 
the vicinity of the track in all layers. The truncatec^ mean value of the energy deposited is then 
divided by the expected track length in the material. 

The specific energy loss by a charged particle depends from the velocity of the particle in 
the material as illustrated by figures [U [2j The relationship is described by the well known 
Bethe-Bloch formula. To obtain an empirical formula for the specific energy loss in the pixel 
detector a fit is made on the dE/dx distributions from low momentum standard model particles 
with well-defined masses (fig. [T]) . 




Figure 1: Measured specific energy loss 
(dE/dx) in the pixel detector as a func- 
tion of momentum for low pr standard 
model particles The bands corre- 
sponds to pions, kaons and protons re- 
spectively. 



Figure 2: Specific energy loss as a func- 
tion of momentum for four different sim- 
ulated R-hadron massesffl 100, 300, 500 
and 700 GeV. 



"This is done to avoid fluctuations from a deposition far out on the Landau tail 



Table 1: Cut flow table: observed data yields and signal expectationM The quoted uncertainties are statistical. 





Data 


Bkg 


Gluino 400 GeV 


Gluino 700 GeV 


Cut level 


# Events 


Cut Eff. 


Total Eff. 


# Events 


Total Eff. 


Total Eff. 


Trigger 


2,413,863 








0.205 ±0.013 


0.219 ±0.009 


Offline EJp^'" 


1,421,497 


0.589 


0.589 




0.200 ±0.013 


0.216 ±0.009 


Primary vtx 


1,368,821 


0.963 


0.567 




0.200 ±0.013 


0.216 ±0.009 


High-pT 


212,464 


0.155 


0.0880 




0.120 ±0.010 


0.129 ±0.007 


Isolation 


32, 188 


0.151 


0.0133 




0.100 ±0.009 


0.105 ± 0.006 


High-p 


21,040 


0.654 


8.7 X 10^3 




0.099 ± 0.009 


0.104 ±0.006 


lonisation 


333 


0.016 


1.4 X 10"'' 


332 


0.067 ±0.008 


0.085 ±0.005 



Equation [T] expresses the most probable value of the specific energy loss for a given velocity 
13. The function is numerically inverted to find the velocity as a function of energy loss and then 
the mass. 

MVV^E/<i.{P) = ^ In (1 ± {p2hY') - PA (1) 

4 Selection 

4-1 Trigger 

It is currently not possible to trigger directly on heavily ionising tracks with the ATLAS de- 
tector. To compensate, the analysis uses a missing energy trigger that relies on the production 
mechanism of the gluinos which also induces associated QCD radiation. A large missing energy 
component from initial state radiation (ISR) jets enables the use of the lowest unprescaled miss- 
ing Et trigger available in the data taking period. The threshold is set to E™'^^ > 70 GeV. The 
overall acceptance in the mass range of interest is around 20%, with a slight mass dependence. 

4-2 Candidate selection 

The selection reduces overall background by requiring offline EJp^'^'^ > 85 GeV, a vertex with at 
least 5 tracks of which one track must have px > 50 GeV and the longitudinal and transverse 
impact parameters below 1.5 mm. High-p-r track candidates must also have three hits in the 
pixel detector, six hits in the silicon tracker and be isolated from other tracks with px > 5 GeV 
and ^Rtrk > 0.25. Furthermore the track must have a momentum p larger than 100 GeV. 
The last cut is based on ionisation. The dE/dx estimate is not entirely flat as a function of 
pseudorapidity, to compensate the track must have a energy loss larger than 1.8 — 0.045|?7| ± 
0.115|r?p - 0.033|r/|3 MeV/g cm . The resulting reduction is illustrated in table [T] for data, 
background and two signal sets. 

4-3 Background estimation 

The background is expected to be entirely due to instrumental misclassification. A full Monte- 
Carlo simulation from first principles would require a very accurate model of the detector feed- 
back to a degree not yet feasible. Instead a data-driven approach is used. By assuming indepen- 
dence between momentum and charge deposition for the background, two independent side-band 
windows {Bkgl and Bkg2 in table l4?3|) can be defined outside the relevant physics region. 

In the Bkgl region momentum is sampled as a function of rj with the requirement that 
dE/dx must be less than 1.8 MeVg~^cm^. In the Bkg2 region the dE/dx spectrum is sampled 
as a function of rj for p < 100 GeV. Tracks in these regions are collected in 2D histograms 
Bkgl{rj,p), Bkg2{r], dE/dx). The two distributions are then sampled randomly 2 million times 





Table 2: Sideband 


regions 




Low dE/dx 


High dE/dx 


Low p 




Bkg2 


High p 


Bkgl 


signal 



for each sampling and the values are combined to form a mass estimate based on the dE/dx 
and the momentum sampled. The resulting distribution is scaled to match the distribution in 
data by matching the sideband region below the signal region m < 140 GeV. The resulting 
background is shown in figure [3l 



5 Results 



After the final selection no significant deviation from the standard model is observed. An upper 
limit is calculated for each mass hypothesis with the CLg method at 95% CL. Comparing with 
expected cross sections produced with Prospino^ a mass exclusion for Split-SUSY gluinos are 
found for gluino masses less than 810 GeV assuming gluino-gluon formation at 10%. Limits for 
all tested masses are shown in figure [H 




Figure 3: Mass distribution for data, background and Figure 4: Upper limit on gluino R-Hadron production 
four simulated gluino R-Hadron signal scenarios scaled cross section^, 

to the Prospino expected cross sections^. 
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